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Abstract 

I examine a possibility to extract the angle 7 of the Cabibbo-Kobayashi- 
Maskawa unitarity triangle by the inclusive B i~ decay. An indepen- 

dent information for the angle is expected from the non-trivial contribution 
induced in uu and cc loops. The contributions induce CP asymmetry which 
has high sensitivity to the angle 7 in the whole dilepton invariant mass region. 
Particularly, in the low dilepton invariant mass region, the sensitivity is also 
recognized in the branching-ratio with any dilepton final states and the lepton 
polarization asymmetry with dimuon final state. In the high dilepton invari- 
ant mass region, the sensitivity is tiny in all measurements with the exception 
of the CP asymmetry, that make them to be good probes to confirm the mea- 
surement of Vtd*Vtb in addition to the present data from — B^ mixing. The 
decay rate and asymmetries are examined in the Standard Model with taking 
into account the long-distance contributions due to vector-mesons as well as 
its momentum dependences that would reduce the long-distance backgrounds 
in the channel. 



*0n leave from P3FT-LIPI, Indonesia. E-mail address : handoko@theo.phys.sci.hiroshima- 
u. ac.jp 



1 Introduction 



The hope that B ^ Xsi~^ i~ decay will be within experimental reach in the near 
future encourage me to consider B i~ decay. Both decays are important 

probes of the effective Hamiltonian governing the flavor-changing neutral current 
(FCNC) transition b ^ q£^i^ {q = s,d) in the Standard Model (SM) as written 
below 

Hes = -4^V,lVtb ICg^'^ [qi.Lb] \iYi\+Cio [q-f^Lb] llYl^i 



-2G 



cff 



qicTfiu-r {R + rhqL) b 



where L/ R= (1 =1= 75)/2, q^ denotes four-momentum of the dilepton, s = q^. Nota- 
tion with hat on the top means it is normalized with mf,. 

Theoretically, the most important interest in the b ^ di~^ i~ decay is, the matrix 
element contains un-negligible terms induced by continuum part of uu and cc loops 
proportional to V*^ Vub and V*^ Vet 0] • These terms should give non-trivial contribu- 
tions in the Wilson coefficient Cg^^ which induce CP violation in the channel. I call 
the contribution as CP violation factor (Cg^^) throughout this paper. I will show 
that it can be utilized to determine the length x and the angle 7 of the Cabibbo- 
Kobayashi-Maskawa (CKM) unitarity triangle in Fig. |1] at once. In this meaning, 
for example, the radiative B 7 decay is not so useful since it gives only infor- 

mation for the length x that have already been measured well in the B^ — B^ mixing. 
On the other hand, in the b — > si~^ i~ decay, Cg^^ is strongly suppressed due to 
the GIM mechanism. Generally, rare B decays are clean processes to extract the 
CKM matrix elements, because non-perturbative effects in the decays are possibly 
tiny, less than few percents as shown in by using heavy-quark effective theory 
approach 0. In the analysis, I also utilize the experiment result of x^ in the B'^ — B^ 
mixing. 

The purpose of this paper is to show a possibility to give an independent mea- 



^ Remark that the approach is reUable only in the low dilepton mass region, but it is sufficient 
to justify the statement since in the high dilepton mass region the perturbative calculation is good. 
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Figure 1: The CKM unitarity triangle on the p — r] plane. 



surement for the angle 7 of CKM unitarity triangle in the SM by observing the 
channel. The calculation is done with taking into account the g— dependence in the 
long-distance (LD) contributions due to the vector mesons 0. It is well known that 
including the g— dependence, which has not been considered in the previous papers, 
will reduce the background due to LD contributions [§]. 

This paper is organized as follows. First I briefly describe the non-trivial con- 
tributions due to resonance and continuum parts of uu and cc loops which induce 
Cg*"^. Next I consider the phenomenology of the contributions and its relation with 
the CKM unitarity triangle. Before going to summary, I analyse the decay rate and 
asymmetries, i.e. forward-backward (FB) asymmetry (^fb), CP asymmetry (^cp) 
and lepton-polarization (LP) asymmetry (.Alp) in the channel. 

2 CP violation factor 

The effective Hamiltonian in Eq. (|l|) describes both inclusive b ^ qi~^ i~ decays by 
replacing q with s or d quark respectively. In the SM, the QCD corrected Wilson 
coefficients enter in the physical decay amplitude above have been calculated up to 
next-to leading order (NLO) for Cg'^^ and leading order (LO) for Cj*^^ [Q, while Cio 
does not receive any correction at all. Some corrections due to continuum parts of 
uu and cc continuums and the resonances of the vector-mesons will enter only in 
the coefficient Cg'^^. Remind that the contribution of uu loop in €7'^^ is suppressed 



due to the GIM mechanism 10 



Before going to give Cg^^, related with uu and cc loops, let me mention the 
operators govern the b —>■ qu^u^ processes, 



01 = {qal^^Lh^) {^'pl^L 

i=l,2 

02 = {qal.Lhp) {^'pl'L 



u, 



j=l,2 



(2) 
(3) 



after doing Fierz transformation. Here, = u, v? = c and the lower suffixes denote 
the color. For q = s, uu loop contribution can be ignored 0, while for q = d the 
situation is quite different. The reason is, both operators above are proportional to 
the CKM matrix element V*,^ V^ih/Vf* Vth if we normalize the amplitude with V^* Vtb 



as usual, then 



0(1) 



u = u 
u^ = c 



for the former channel, and 



0(1) , u' = u 
0(1) , u' = c 



(4) 



(5) 



for the later one. A is a parameter in Wolfenstein parametrization of CKM matrix 
1^ and the world average is A ~ 0.22 

Involving the continuum as well as resonance parts of uu and cc loops and NLO 
QCD correction into the calculation gives 



ri eS _ NLO 
L'9 — <-^9 



1 + ^co(s) 



(6) 



where 



09^°'^(S) 



V* T/ 

g[mc, s) - -———g{mu, s) 



VtlVtbJ''' ' V,\Vib 
X (30i + O2 + 3O3 + O4 + 3O5 + Oe) 

-\g{\,s) (403 + 404 + 305 + 06) 

-i(?(0, s) (O3 + 3O4) + ^ (3O3 + O4 + 3O5 + Oe 

(3O1+O2+3O3 + O4 + 3O5 + O6) 



(7) 
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X 



^4 ^tfe / y=v, 



Vtg Vtb v=p,uj 



The readers should refer |E| for Cg^^'^ and 



2s(l + s)(l -2s) , . 5 + 9s-6s2 



(9) 



3(1 -s)2(l + 2s) 6(l-s)(l + 2s 

represents the 0{as) correction from the one gluon exchange in the matrix element 
of Og. The function g{m^i,s) which describes the continuum part of pair 
contribution is 

16 2 



-°lnp^ 
9 \u 



- ln(m,,i) H , 

9 ^ ''^ 27 9 s 



9l^ + ^ 



1 - 



e 1- 



4m 



l + Vl-4<./S , 
m , — z TT 



1 



1 - 4m^,/'5 



+0 f ^H^L _ 1^ 2 arctan 



^7(0,S) 



s 



4m2,/s - 1 



(10) 



In , 

27 9 \ii 



4, . 4 . 
- m s H — z vr 
9 9 



(11) 



In the resonance part Cg^^^, I put the relative phase to be zero because of unitarity 
constraint in the Argand plot of the transition amphtude [0. -FV(s) is the Breit- 
Wigner resonance form 



Fv s) 



s — rhy + i rhv Ty 



(12) 



fv{s) describes the momentum dependence of coupling strength of vector interaction 
in 7 — transition, i.e. 







and has been derived as follows 

fv{s) 



MO) Pv^ ^ 



(13) 



(14) 



under an assumption that the vector- mesons are bound-states of the pair Here, 

,2 



m: 



47r2 s 



5 s s 
-4 - — ^ + 41 + 



3 m, i 



2m^ 



— 1 arctan 



4m2^/s - 1 



(15) 



Parameter 






A/pftor- 











r 


OJ 


lb 

T 


ih' 

T 


tb" 


tb'" 


My (MeV) 


768.5 


781.94 


3096.88 


3686.00 


3769.9 


4040 


Te+e- (keV) 


6.77 


0.6 


5.26 


2.14 


0.26 


0.75 


Tv (MeV) 


150.7 


8.43 


0.087 


0.277 


23.6 


52 


-0.00243 


-0.00243 


-0.02734 


-0.01463 


-0.01374 


-0.01153 


Pv 


0.01339 


0.01387 


0.01577 


0.01830 


0.01885 


0.02080 


/v(0) 


0.00662 


0.00202 


0.01795 


0.01487 


0.00536 


0.01010 



Table 1: The experimental (central) values for each vector- meson under considera- 
tion (upper table) and the determined constants under these values (lower table). 



is obtained from a dispersion relation involving the imaginary part of quark-loop 
diagram, while Py and Py are the subtraction constants. Kinematically the above 
interpolation equation of fy is valid only for < s < rhy region. For s > rhy 
region I take same assumption with P], that is fy{s > Thy) = fv{my). In principle, 



the ratio in Eq. (jM]) should be obtained from the known data of V production 
cross-section by off-shell and on-shell photons. 

The subtraction constants in Eq. (|l^) are written in the lower table of Tab. (|l]) 
for each vector meson. The results are determined by using the data in the upper 
table, putting my ~ (2m„i) and Py = 0.043 for all Vs. Unfortunately, there is no 
data of photoproduction for higher excited states of ip, so let me use same average 
value |/y(0)//y(m^)| ~ 0.35 for V = , , ^j" , ijj'" and |/y(0)//y(m^)| ~ 0.92 for 
V = p,uj which fit the data on p,uj and ip 0. This fact is also the reason why 
other resonances higher than ip"' are not considered here. Otherwise, fy^rhy) can 



be obtained from the data on leptonic width |11], that is 



/v(^v) = §^f(i^-^^n, (16) 

16 vra^ 



then /y(0) would follow respectively as written in Tab. |l[ 

From Eqs. (|^) and (§), it is obvious that V*g Vub/V^* Vtb would induce CP violation 
in the channel. For convenience, these terms can be collected as {y*^ Vub/V^* Vtb^ Cg'"^{s), 
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Parameter Value 



mw 80.26 ± 0.16 (GeV) 

mz 91.19 ±0.002 (GeV) 

rriu 0.005 (GeV) 

ma 0.139 (GeV) 

1.4 (GeV) 

nif, 4.8 (GeV) 

mt 175 ± 9 (GeV) 

me 0.511 (MeV) 

105.66 (MeV) 

nir 1777l[J:^? (MeV) 

II btli (GeV) 

Ag^D 0.214;°:°^^ (GeV) 

c^QED 129 

«,(mz) 0.117 ±0.005 

sin^^^ 0.2325 

niBO 5279.2 ± 1.8 (MeV) 

tbo 1.28 ± 0.06 (ps) 

rjQCD 0.55 



fB^BB, 173 ±40 (MeV) 
BiB^XJu) (10.4 ±0.4)% 
Xd 0.73 ±0.05 

Table 2: The values of parameters used throughout the paper. 



with 

C9^P(S) = (3Ci±C2±3C3±C4±3C5±C6 



X 



16 TT^ 

g{mc, s) - g{rhu, s) ^| J2 ^vis) - J2 ^v(^) 



(17) 



Now let me derive some relations in the CKM unitarity triangle and give numer- 
ical calculation for the auxiliary functions defined above. Especially it is worthwhile 
to see how large the contribution of Cg*"^. Using Wolfenstein parametrization 0, 
one can rewrite the CKM factor as 



-^-\ l±>-2rcos7 ' ^ = ^ • (18) 



(e-'^ - r) 

[ AVe-*^ , q = s 

As mentioned before, from Eqs. (^) and ([ISD it is obvious that in the b si^ i~ 
decay, uu loop and then Cg^^ is less important and negligible. In general one 



Vu*d ^ub ^ CP 




Figure 2: The magnitude of Cg^^ with (sohd thick curve) and without (sohd thin 
curve) resonances for (7,x) = (48.6°, 0.778). The upper, middle and lower dashed 
curves show the magnitude without resonances for same x and 7 = 0°, —90°, 90°. 

must treat r and 7 as free parameters, while x must be determined by the data of 
— mixing. However, in the SM the unitarity triangle is satisfied in a good 
aproximation, so one can relate r and x each other as below 



x^ — sin 7I + cos 7 . 



(19) 



Here x is determined by the experimental value of Xd in the B^ — B^ mixing, that is 



X 



with 12 



GfV (6 7r2) TUBa ^w'^ Tb^ VQCD Jb/ Bb^ |-^AB=2| 

4 — 15 + 12 Xt^ — — 6 Xt^ In Xt 



1/2 



(20) 



F/^B= 



and Xt = {mt/mw)^. For m^^ and r^^, I use m^o and r^o. 

The magnitude of (KTcZ ^fe/^d ^fc) C'g*"^ ^ function of s is depicted in Fig 
Here I use the following values for the Wilson coefficients 

Ci = -0.2404, C2 = 1.1032, C3 = 0.0107, C4 = -0.0249, C5 = 0.0072, 



(21) 
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Cq = -0.03024, Cj"^ = -0.3109, Cg = -0.1478, Cg^^^ = 4.1990, do = -4.5399, 

that are obtained by using the central values in Tab. (Q), The solid curves show 
the magnitude for (7,x) = (48.6°, 0.778) that is the best fit in the SM up to now 
and equivalent to (p, r/) = (0.3,0.34). The upper, middle and lower dashed curves 
show the magnitude without LD effects with 7 = 0°, —90°, 90°. It is obvious that 
the contribution is significant, about ~ 20% of Cg^^*^ at low s region (s < 0.4). 
Now, I am ready to analyse the decay rate and asymmetries in the channel. 

3 Decay rate and asymmetries 



The double differential decay rate for semi-leptonic B Xq 
the lepton and light quark masses, is expressed as 



decay, involving 



dsdz 



/~i eff 



\c 



10 1 



mo 



s + 771': 



+ 
+4 



a 



eflf 



IC- 



10 



[1 - rhlY ~ - u{s)^ 1- 



6 



2 1 + 2m^/s 



X 



1 — — 171^ + — s (8 + s + s) + u{sf' (1 + m^) 



-8RefCg'=f^V Cj""^ 



1 + 



2 m? 



s{l + ml)-{l-rnlf 



+4Cio [Re (Cg^^)* s + 2^7^^ [l + rn])] u{s) z] 



(22) 



where u{s) = — (1 + rfig)^] [s — (1 — m^)^], z = cos 6 is the angle of mea- 
sured with respect to the 6— quark direction in the dilepton CM system and the 
normalization factor, 

2 



Bo = B{B ^ XJu 



3 a' 



1 



(23) 



is to reduce the uncertainty due to 6— quark mass. In the preceding notation, the 

td ^tb\ /l^cbl ~ A^a;^. fifhc) is the phase space function 



CKM factor would read I V;:i VlK' - ■ 



for V{B Xcdv) in parton model, while nijhc) accounts the 0{as) QCD correction 
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dB 

ds 




Figure 3: Left : differential BR for e~ with (thick curve) and without (thin curve) 
resonances for (7, x) = (48.6°, 0.778). The upper, middle and lower dashed curves 
show the SD contribution for same x and 7 = 0°, 90°, —90°. Right : same with the 
left for 7 = 48.6° and x = 1.368 (solid curve) and 0.631 (dashed curve). 



to the decay. Writing both functions explicitly. 



1 — 8 m\ + 8 m\ — — 24 In m, 
2as{mh) 



37T 



3 / 2 31' . 



(24) 
(25) 



and using the values in Tab. 0, f{rhc) = 0.542 and K{rh^ = 0.885. 



3.1 Decay rate 

Now, let me consider the dilepton invariant mass distribution of the differential 
branching-ratio (BR). It can be obtained by integrating Eq. (^) over the whole 
region of variable z, 

(26) 



dB{s) _ ^ d^B{s,z) 
ds J-i dsdz 



This gives 
dB{s) 



ds 



u(s) < 6 



efF 



- \c- 



10 



1 — s + mt 



+ 



cff 

Og 



+ \c- 



10 1 



;i - m^)2 + s (1 + mj) - 2 + u{sf 



2 m? 



10 



+4 



2 1 + 2m|/s 



2(l + mJ)(l-mJ)2- (l + 14m^ + mj)s- (l+mj)s' 
2mT 



+ 12Re(C9^'^)* C, 



off 



1 + 



'1 



2\2 



(27) 



The distribution of differential BR on dilepton invariant mass for B —>■ e+ e~ 
is given in Fig. ^for various values of 7 and x. High sensitivity on x is mostly coming 
from the CKM factor in Eq. (|23|), while the sensitivity on 7 seems significant only in 
the low s region. Then, in the high s region the differential BR may be a good test 
for X and makes good the loss of — B^ mixing due to the theoretical uncertainties 



B). Anyway, I have checked 



in the treatment of hadron matrix element (^B O^O 
that the momentum dependence of resonances in the channel is not as large as the 
result in There is only ~ 4% reducement compared with using fvi^v) for all 
region. This is because the 1/s suppression in Eq. (p^ [Q. 

Next, I am going on examining some measurements that are sensitive on 7 and 
less sensitive on x. This may be achieved by considering the asymmetries and 
normalizing them with the differential BR to eliminate the CKM factor. 

3.2 Forward- backward asymmetry 

First, I provide the FB asymmetry. The normalized FB asymmetry is defined as 
follows M] 



Afb{s) 







dsdz 



^^ d^B{s,z) 

^1 ^ dsdz d^FB(s)/ds 



1 d^B{s,z) ^ /-o d^B{s,z) dB{s)/ds 
dsdz J-i dsdz 



(28) 



Then, after integrating Eq. ( jn\ ) properly, the nominator reads 
d^FB(s) 



-4S„Wl 



4 mi 



uisYC 



10 



Re (Cg"^)* 5 + 2^7"^ (1 + 



(29) 



ds V s 

In Fig. ^, I plot the FB asymmetry for B e+ e~ with and without the 

resonances in the left figure, while in the right one with varying 7 and keeping x to 
be constant. 
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Figure 4: Left : FB asymmetry for e+ e~ with (thick curve) and without (thin curve) 
resonances for (7,0;) = (48.6", 0.778). Right : same with the left for same x and 
7 = 90° (sohd hne) and 0° (dashed hne). 

3.3 CP asymmetry 

Doing same treatment as in the amphtude level, the normalized CP asymmetry 
can be written simply as 

di3(s)/ds - dg(s)/ds -2 d^cp(s)/ds 



^cp(s) 



(30) 



dB{s)/ds + dB{s)/ds dB{s)/ds + 2 d^cp(s)/ds 

where B and B denote the BR ofb^qi^£~ and its complex conjugate b — ^ g£+ i~ 
respectively. For convenience, Cg'^^ is divided to be two terms according to the factor 
Vuq Vub/V^l Vtb as below 



T/* T/ 

n _ I ^uq ^ub „ CP 

Then, the result for the differential CP asymmetry is 
d^cp(s) 4„ 4rfi^ 



(31) 



ds 



■u(s) Im 



X < Im 



+6 Im 



/=y* CP 



:i - fnlf + s(l + m1)-2 + u{s) 



,2 m? 



+6m^ (1 — s + ?fig) 



^ eff CP 

07 09 



1 + 



2 m} 



{l-rnlf-{l + rnl)s . (32) 
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CP 
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Figure 5: Left : CP asymmetry for e~ with (solid thick curve) and without (sohd 
thin curve) resonances for (7,x) = (48.6°, 0.778). The upper and lower dashed 
curves show the SD contribution for same x and 7 = 90°, —90°. Right : X as a 
function of 7 where ^cp = 0. The dashed lines denote the upper and lower bounds 
of X in the SM. 

In Fig. 1^, I give the distribution of ^cp (left figure) with (solid thick line) and 
without (solid thin line) resonances. It is easily understood that in the present case, 
the dependence on 7 is large, because of the appearence of factor rsin7 from the 
CKM factor in Eq. (|32D . Moreover, it is clear that ^cp will be non-zero if the 
imaginary part of Eq. (pIS]) is non-zero. Anyway, a condition that Acp = for q = d 
in the SM is satisfied by the following equation. 



2 • 2 

X = sm 7 



1 + 



3 + 4 cot 7 



(33) 



by using Eq. ([T9|). The right figure in Fig. |] is plotted based on this equation. As de- 
picted in the figure, there are still allowed regions of 7 where ^cp = 0. Notice again 
that from Eq. (0), Acp{B ^ X, £+ £-) would be ~ 5% of AcAB ^ i+ £") due 
to the suppression of A^. 

3.4 Lepton-polarization asymmetry 

Until now, all of the measurements are less sensitive to the lepton mass. Next, 
I provide the LP asymmetry which must be considered for heavy dilepton final 
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Figure 6: Left : LP asymmetry for r"*" t" with (solid thick curve) and without 
(sohd thin curve) resonances for (7,x) = (48.6°, 0.778). Right : same with the 
left but for The dashed curves show the SD contribution for same x and 

7 = o°,90°,-90°. 



state and is proposed first in [|T4| for B ^ X^t^ t . Generally, the normalized LP 
asymmetry is given as 

dB{s, n)/ds - dB{s, -n)/ds dAp(s)/ds 



Alp(s) 



(34) 



dB{rs,n)/ds + dB{rs,-n)/drs dB{s)/ds 
with n is a unit vector of any given spin direction of i~ in its rest frame. Then, for 
the longitudinal polarization, that is n has same direction with the momentum of 

dAMs) 



ds 



+Re(C9'=«f)' 



MS 



m 



2\2 



+ s(l + mg^) -2s^ } 



(35) 



The distribution is depicted in Fig. As shown in the figure, for r+r^ final 
state the sensitivity on 7 is tiny. The reason is, for r+r^ final state the distribution 
starts appearing in the region higher than s = (4m^), on the other hand, generally 
the high sensitivity on 7 is expected in the low s region. So it may be interesting to 



consider the transversal polarization that has different structure [|14[ . Unfortunately, 
I have checked that the transversal lepton polarization asymmetry is too small for 



0.2 0.4 0.6 0.8 1.0 

s 
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light dilepton like /x"*"// , then one must consider t'^t again that the distribution 
is limited for higher s region. 

4 Summary 

I have shown how to extract the angle 7 of CKM unitarity triangle by the inclusive 
B — > Xd i~ decay and the data of in the SM. As the results, finally 1 can make 
some points as below. 

1. For low s region, i.e. 0.1 < s < 0.3, 5 ~ 10% discrepancies in B{B — > X^ie'^ e~) 
and Ahp{B Xd/J^'^ IJ>~) may be good signals for the CP violation factor 
defined here. 

2. According to the fact that the sensitivity on 7 in the high s region, i.e. s > 0.6, 
is tiny, a measurement of x may be done by exploring one of the measurements 
discussed in the present paper in addition to the present data of Xd in — B^ 
mixing. 

3. CP violation asymmetry in the channel should measure the dependence on 7. 
It will, at least, be a good probe to determine the sign of the angle 7. 

To conclude, the measurements of the decay rate and asymmetries in B ^ 
Xdi^ i~ decay will provide an independent information for 7. The information 
is a crucial test of CKM unitarity as well as leading to the discovery of unitarity 
violation. 
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